Mott-insulator states of ultracold atoms in optical resonators 
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We study the low temperature physics of an ultracold atomic gas in the potential formed inside 
a pumped optical resonator. Here, the height of the cavity potential, and hence the quantum state 
of the gas, depends not only on the pump parameters, but also on the atomic density through 
a dynamical a.c.-Stark shift of the cavity resonance. We derive the Bose-Hubbard model in one 
dimension, and use the strong coupling expansion to determine the parameter regime in which the 
system is in the Mott-insulator state. We predict the existence of overlapping, competing Mott 
states, and bistable behavior in the vicinity of the shifted cavity resonance, controlled by the pump 
parameters. Outside these parameter regions, the state of the system is in most cases superfiuid. 
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Ultracold atomic gases in optical lattices offer the un- 
precedented and unique possibility to study paradigmatic 
systems of quantum many-body physics [3, 0]. These 
systems allow one to realize various versions of Hubbard 
models Q, a prominent example of which is the Bose- 
Hubbard model 0] , exhibiting the superfiuid (SF) - Mott 
insulator (MI) quantum phase transition [f|. The real- 
ization of the Bose-Hubbard model with ultracold atoms 
has been proposed in the seminal theoretical work in 
Ref. and demonstrated in the milestone experiments 
in Ref. ■ Several aspects and modifications of the SF - 
MI quantum phase transition (or crossover 0]) are object 
of intense studies [2| • 

Optical lattices in free space are not affected by the 
presence of the atoms. This scenario is, however, strongly 
modified when the atoms move in the optical potential 
which is formed inside a pumped resonator: Here, the 
atoms interact with the cavity mode while the cavity 
field, determining the optical lattice, may critically de- 
pend on the density of the atoms [9|, [lfj. Several re- 
cent studies address Cavity Quantum Electrodynamics 
(CQED) with cold atoms. CQED techniques were used 
to measure pair correlations in the atom laser [Til ] . and 
have been proposed for characterizing quantum states 
of ultracold matter (l2| . Self-organization of atoms in 
transversally pumped cavities was observed in and 
theoretically described in [lj] . Bragg scattering of atomic 
structures inside optical resonators has been investigated 
in 15|. Most recently, Bose-Einstein condensed atoms 
have been loaded inside cavities [lrjj ]. This experimen- 
tal progress calls for theoretical development of CQED 
combined with many-body physics. 

In this Letter we determine the ground state of ultra- 
cold atomic gases in the optical lattice of a cavity. The 
cavity is driven by a laser, and the atoms shift the cavity 
resonance, thus affecting the intracavity field amplitude, 
which in turn determines the depth of the cavity poten- 



tial and the ground state of the atomic gas itself. The 
problem is hence highly non-linear, as the optical lat- 
tice and state of the atoms have to be evaluated in a 
self-consistent way. The derivation of the corresponding 
Bose-Hubbard model for few atoms has been discussed 
by Maschler and Ritsch in Refs. 17]. In this Letter, we 



derive the Bose-Hubbard model in an appropriate ther- 
modynamic limit. We study its ground state applying the 
strong coupling expansion [l8| to calculate the bound- 
aries of the MI states, determined by the dependence on 
the parameters of the system: pump strength and fre- 
quency, density of atoms, and chemical potential. 

Our model consists of bosonic atoms confined in a ID 
trap inside an optical resonator of a fixed length driven 
by a laser field. The atomic dipole transition is far-off 
resonance from the cavity mode, which induces a di pol e 
potential acting on the atoms. Using the notation of 
the single-particle Hamiltonian reads: 



17j, 



Ho = Vh [Uq cos 2 (fci) — A c ] fiph — ihr] (d — d^) . (1) 



Here, p, x and m are the atomic momentum, position, 
and mass, rj is the amplitude of the pump at frequency 
u)p, A a = uj p — uj a and A c = lu p — oj c are the detunings of 
the pump from atom and cavity frequencies, k — lu c /c is 
the mode wave vector, d^ and d are the creation and an- 
nihilation operators of a cavity photon of energy Huj c , and 
fipb = d^d is the number of photons. The depth of the 
single-photon dipole potential is Uq — <7g/A Q , where go is 
the atom-cavity mode coupling. The many-body Hamil- 
tonian is obtained from Eq. {1} including the atomic con- 
tact interactions; it is conveniently represented in second- 
quantized form with the atomic field operators ^(x), 
^f'(x) obeying the bosonic commutation relations. We 
assume the bad-cavity limit, where the resonator field 
reaches the stationary state on a faster time scale than 
the one of the atomic dynamics, and eliminate the cav- 
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ity field from the equations of the atomic operators. In 
this limit the amplitude of the intracavity field depends 
non-linearly on the atomic fields through the operator 
y = J dx cos 2 (kx)W (x)$> (x) , and reads 



a(j>) = 



n 



k - i(A c - UoS>) 



(2) 



where k is the cavity damping rate. Correspondingly, the 
Heisenberg equation for the atomic field operator reads 



: {^(x),n ]-iC(y,x), 



(3) 



where Ho = 



fdx&(x) (- 



with u being the strength of the contact interaction, and 

60, x) = U cos 2 (kx)a"<(y)$(x)a(y) (4) 

which arises from keeping track of the correct normal or- 
dering of atomic and photonic field operators. Starting 
from Eq. |3j , the derivation of the Bose- Hubbard Hamil- 
tonian is not straightforward due to the form of opera- 
tor p)). This is evident when expanding the atomic field 
operators, assuming the validity of the tight-binding ap- 
proximation (TBA) and the occupation of the lowest en- 
ergy band: $f(x) — J^L w(x—Xi)bi, where hi and w(x— Xi) 
are the atomic annihilation operator and Wannier func- 
tion at site i, respectively. The Wannier functions depend 
on y and, therefore, on the number of atoms N or, equiv- 
alently, on the atomic density, which in turn depends on 
Wannier functions. Moreover, the commutation relation 
[bi, bj] = Sij is valid only in the lowest order in the expan- 
sion in 1/N. In effect, the Wannier expansion must be 
performed self-consistently in the thermodynamic limit, 
by letting TV and the cavity volume to infinity, keeping 
finite the number of atoms per potential site. Addition- 
ally, we impose Uq — uq/N and r\ = rjo^/N, where uq 
and r/o are constants, which corresponds to keeping the 
depth of the cavity potential V — hUon p h constant as 
N is increased. The Bose-Hubbard Hamiltonian H is 
obtained discarding couplings beyond nearest-neighbor. 

Its rescaled form H = H/U, with U the strength of the 
on-site interaction, reads 



H 



-tB 



1 



is expressed in terms of ( — A c — uoJofio and of the 
coefficients U — u/2 J dx \w(x)\ 4 , Eg — f dxw(x — 
xi)(—h 2 /2m)(d 2 /dx 2 )w(x — xi + i) and Jg = J dxw(x — 
xi) cos 2 (kx)w(x — xi + g), with I = 0, 1. Note that N, and 
hence the atomic density in the homogeneous case, is a 
conserved quantity since [N, H] = 0. In deriving Eq. ([5|) 
we have used that J\ <C Jo- The higher order terms in 
J\B, describing long-range interactions, have hence been 
neglected. We note that the parameters jx and t depend 
on the atomic density through the Wannier functions, 
and at the same time determine the state of the system, 
and in particular the density: This is a genuine CQED 
effect, where the non-linearity of the coupling between 
photons and atoms depends on the atom number. As 
a consequence, the atomic density in this system is not 
determined by the chemical potential alone. 

From Eq. (|5|) we determine the parameter regimes of 
the MI states using the strong coupling expansion [l8| . 
Here, the boundaries of the MI regions are determined 
by comparing the energy of the MI state, given by no 
atoms at each site of the periodic potential, with the cor- 
responding energies of the excited states with one ad- 
ditional or missing particle (particle and hole states). 
This procedure involves the evaluation of the Wannier 
functions for all cases. As the coefficients of the Bose- 
Hubbard Hamiltonian both depend on and determine 
the atomic density, the Wannier functions have to be 
calculated self-consistently. This is done by solving the 
non-linear equation in presence of the potential V — 



hUonp h = hu ri 2 / 



where Jq is an 



1) 



fiN, 



(5) 



+ (A c - u J n ) 
integral of Wannier functions. This equation cannot be 
solved by iteration, as one encounters periodic doubling 
bifurcations and deterministic chaos. We solve it numer- 
ically by checking for self-consistent solutions, using the 
Gaussian approximation of the Wannier functions, and 
thus approximating w(x) ~ exp(— x 2 /2a 2 )/ (y/na) 1 ^ 2 
where a is the parameter to be determined [19( . In terms 
of the dimensionsless quantity y = k 2 a 2 , giving the ex- 
tension of the Gaussian wave packet in units of the cavity 
mode wavelength, the problem can be reduced to solv- 
ing self-consistently the equation Jo{y) = | (1 — d a e~ v ), 
where d a is the sign of the detuning A a . For a given set 
of parameters, multiple (bistable |20() solutions appear 
when the number of photons is maximum, namely when 
the denominator of Eq. ([2|) is minimum, which occurs at 
the shifted resonance 



where N = J2i ™i — J2i K^i is the atom number operator 
and B — b\hi-i + h.c. is the hopping term. The term 

ft = + contains the rescaled chemical potential, 
where the second term is a constant in the thermody- 
namic limit. The tunneling parameter 



E 1 KrfUpJx 

W ~ u {n 2 + c 2 ) 



(6) 



A c - uo Jmo = 0. 



(7) 



Since the sign of uo is determined by the detuning A a , 
Eq. ([7]) allows for real solutions only when A Q and A c 
have the same sign: Then, the resonance condition de- 
pends on the number of atoms. Correspondingly, the cav- 
ity is driven at resonance, the number of photons reaches 
the maximum value n p h = rf 1 /k 2 , and the cavity poten- 
tial is the deepest. An important distinction must be 
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made between the cases A a > (Uq > 0) and A a < 
(Uq < 0): In the first case, the potential minima are at 
the nodes of the standing wave, where the intracavity 
field vanishes. Strong localization of the atoms at these 
points implies that the coupling of the atoms with the 
field is minimum, Jo — > 0. The quantum fluctuations 
give rise to a finite coupling, determining the quantum 
state. On the contrary, when A < the potential min- 
ima are at the antinodes of the standing wave, where the 
intracavity field is maximum. Strong localization of the 
atoms at these points implies strong coupling with the 
field, with Jo — ► 1. In this regime, CQED effects are 
expected to play a dominant role. 

We plot the boundaries of the resulting Mott states in 
the fx — i]^ 1 plane, i.e., the effective chemical potential 
and the inverse of the pump strength. Here, large pump 
strengths correspond to deep optical potentials, hence 
to vanishing tunneling, t — * 0. The physical system we 
consider is a gas of 87 Rb atoms with scattering length 
a s = 5.77 nm, whose dipolc transition at wavelength A = 
830 nm couples to the mode of a resonator at decay rate 
k = 2tt x 100 kHz. The potential has transverse size 
A y = A z = 30 nm and K sites in the longitudinal axis. 
We evaluate the "phase diagrams" for K = 50 — 10000 
at fixed number of atoms N, scaling N so to keep the 
atomic density constant. The results for the Mott zones 
agree over the whole range of values, so in the figures we 
report the ones obtained for K = 50. 

We first discuss the case in which the detunings A a 
and A c have different signs, i.e., far from bistability (2pj 
when Eq. ([7|) is not fulfilled. As expected, there is a 
peculiar difference between the cases A a > (atoms at 
the nodes) and A a < (atoms at the antinodes): When 
A a < 0, in the tight- binding regime Jo — > 1 and thus 
V oc 1/(k 2 + (A c + n u ) 2 ). Hence, the dependence on 
the atomic density is strongest. In Fig.QJa) the MI zones 
are displayed. The shapes are similar to the standard 
ones [18| . apart from the region n — > 00 where, despite 
being out of the bistable regime [2fJ, the lobes consider- 
ably overlap. This effect is due to the competition be- 
tween the non-linear coupling to the cavity field, giving 
the depth of the potential, and the strength of the on- 
site interactions, affecting the number of atoms per site 
no- In the other case, A a > 0, one has Jo <C 1. For 
77 — > 00, then Jo — > 0, the cavity potential depth is al- 
most independent of no and one obtains the standard 
Bose-Hubbard Hamiltonian. For large but finite values 
of 77, however, Jo is finite and the dependence of the co- 
efficients on no becomes relevant. Figure QJb) shows the 
"phase diagram" for \Uo\ = 2k and A c = 0. Here, the MI 
regions exhibit a regular behavior at 77 — * 00. As 77 is de- 
creased they start to overlap and become disconnected. 
This behaviour introduces two new critical points at the 
tips of the disconnected regions, whose nature will be 
studied in future works. We note that the MI zones en- 
ter the region of negative \x. The minima of \x are at the 



pump values where the on-site repulsion is balanced by 
the effective potential V. 
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Figure 1: Boundaries of MI regions as a function of the 
rescaled chemical potential jl and the inverse of the pump 
strength 77 (in units of k) for (a) A c = 2k and Uo = —2k and 
(b) A c = and Uo = 2k; no denotes the site occupation in 
the ID cavity lattice potential of K = 50 sites. The dashed 
lines show the boundaries of the zones which are hidden. 

We now consider the situation when A Q and A c have 
the same sign, such that Eq. (J7J) may have real solutions. 
From Eq. (J7JI, for n = 1 and Uq = —k we find bistable 
behavior for J = |A c |/50k and J sufficiently close to 
1, which is fulfilled for instance for A c = —45k. The 
corresponding diagram is displayed in Fig. [2] The in- 
set shows the potential V as a function of 77 for no = 1. 
Here, one encounters a bistability point while lowering 
the pump intensity, where the cavity field potential dis- 
continuously jumps to a second branch with |V| <SC E r , 
E r = h 2 k 2 /2m being the recoil energy. The first branch 
corresponds to the left Mi-region of the phase diagram 
at no = 1. In the second branch, instead, the TBA is not 
valid, hence most probably the atomic gas will no longer 
be in the lowest band of the cavity potential, and rather 
definitely no longer in a MI state. Using both Gaussian 
and Wannier functions we have verified that our treat- 
ment breaks down as soon as the system goes out from 
the MI region on the left of Fig. [2l This instability has 
the character of a first order transition. The right MI re- 
gion in Fig. [3 is found by applying the theory of [18J. It 
occurs at values of 77 for which \ V\ is in the second branch, 
and is thus of dubious validity, since here the TBA breaks 
down. Instability leads here apparently to population of 
higher Bloch bands; most probably the true ground state 
in this regime is SF (BEC in a very weak lattice poten- 
tial). Let us conjecture that we could find parameter 
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regions where the cavity potential for both branches of 
solutions would support the TBA. In that case, at a given 
density there would exist two stable values of the tunnel- 
ing and onsite-interaction matrix elements. Then, for a 
fixed ft we would have two possible phases, of which only 
one will be energetically favorable, but both being by 
construction stable with respect to small perturbations, 
such as single particle or hole excitations. 




Figure 2: Phase diagram showing the MI regions for A c = 
—45k and Uo = —re. Inset: V (in units of E r ) as a function 
of 7] for no — 1. The bistability causes the system to go out 
of the left MI region with no = 1. The right Mi-region for 
no — 1 is an artifact of the theory [l8l [. as here the TBA is 
invalid. 

So far we have considered a fixed number of atoms. If 
instead the system is coupled to an atomic reservoir, and 
the number of atoms is hence not fixed, then a change 
of the system parameters can lead either to continuous 
or to abrupt changes of the atomic density, and hence 
of the ground state, analogue to second- and first-order 
phase transitions, respectively. When the parameters are 
such, that the system does not exhibit bistability, then ft, 
re/77 an d the atomic density determine uniquely the coef- 
ficients in the Bose- Hubbard model and the ground state. 
In the case of overlapping MI zones, as in Fig. [Ija), the 
system would relax to the state with the density which 
minimizes the energy. Nonetheless, by slowly changing ft, 
corresponding to spanning the phase diagram along the 
ordinate, we expect hysteresis in the atomic on-site den- 
sity, since the states inside the Mott correspond to local 
energy minima. Note that such local changes of on-site 
ft are also encountered when the system is inhomoge- 
neous, for instance, in presence of a harmonic trap. The 
atomic density will exhibit then the characteristic "wed- 
ding cake" form (cf. @, @, Q). Sufficiently slow changes 
of the trap frequency should then lead to hysteresis in the 
"wedding-cake" shape. Applications of this effect could 
include many-body quantum switches and generation of 
coherent superpositions of Mott states. 

This letter refers to the case, in which the atomic den- 
sity globally affects the cavity field. Situations, when 
the atoms may affect locally the cavity field, can be 
found in multimode resonators 2l|, |22j. In these sce- 
narios one could find features typical of phonon-likc 
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